Molecular superfluid phase in one-dimensional multicomponent fermionic cold atoms 
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We study a simple model of iV-component fermions with contact interactions which describes fermionic 
atoms with N = 2F + 1 hyperfine states loaded into a one-dimensional optical lattice. We show by means 
of analytical and numerical approaches that, for attractive interaction, a quasi-long-range molecular superfluid 
phase emerges at low density. In such a phase, the pairing instability is strongly suppressed and the leading 
instability is formed from bound-states made of N fermions. At small density, the molecular superfluid phase is 
generic and exists for a wide range of attractive contact interactions without an SU(iV) symmetry between the 
hyperfine states. 
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Because of rapid progress in recent years, cold atom sys- 
tems have become a major field of research for investigating 
the physics of strong correlations in a widely tunable range 
and in unprecedentedly clean systems Jl]]. Ultracold atomic 
systems also offer direct access to the study of spin degen- 
eracy since the hyperfine spin F can be larger than 1/2, re- 
sulting in 2F + 1 hyperfine states. In non-magnetic traps, 
such as optical traps, this high degeneracy might give rise to 
novel exotic quantum phases. The superfluid state of opti- 
cally trapped alkali fermions with hyperfine spin F > 1/2 has 
been studied with an emphasis on the general structure of the 
large-spin Cooper pairs |2j]. The spin-degeneracy in fermionic 
atoms is also expected to give rise to more complex superfluid 
phases. In particular, a molecular superfluid (MS) phase might 
be stabilized where more than two fermions form a bound 
state. Such a non-trivial superfluid behavior has already been 
found in different contexts. In nuclear physics, a four-particle 
condensate — the a particle — is favored over deuteron conden- 
sation at low densities J3| and it may have implications for 
light nuclei and asymmetric matter in nuclear stars 0]. This 
quartet condensation can also occur in semiconductors with 
the formation of biexcitons A quartetting phase, which 
stems from the pairing of Cooper pairs, has also been found in 
a model of one-dimensional (ID) Josephson junctions |@]. A 
similar phase has also been reported in exact-diagonalization 
calculations of the two-dimensional t-J model at low doping 
0]. More recently, the emergence of quartets and triplets 
(three-fermion bound states) has been proposed to occur in 
the context of ultracold fermionic atoms (H lU [Tol IT ] . 



In view of this increasing interest in the formation of com- 
plex superfluid condensates, it would be highly desirable to 
have at one's disposal a simple paradigmatic A-component 
fermionic model which displays this exotic physics. It will be 
shown in this letter that such a model is provided by the ID 
A-component fermionic Hubbard model with attractive con- 



tact interaction: 

H = -t5^[4 )i Ca,i+i + H.c] 



u 



(1) 



where cl , is the fermion creation operator corresponding 



„t 

to the A = 2F + 1 hyperfine states a = 1,...,N and 
m = J2 a c a i c a,i i s tne density at site i. Model ([TJ dis- 
plays an extended U(A)=U(l)xSU(A) symmetry and it has 
been recently introduced in the context of ultracold fermionic 
atoms [ 12]. A possible experimental realization of this model 
([T]i for A = 3 would be a system of 6 Li atoms loaded into 
a ID optical lattice with a carefully tuned combination of ex- 
ternal magnetic and optical fields to make three internal states 
exhibit SU(3) symmetry 11311 . The A = 4 case might also 
be relevant to the optical trap of four hyperfine states of 40 K 
(F = 9/2 atoms) H. The SU(A) symmetry of Eq. © has 
an important consequence since, when A > 2, even for U < 
there can be no pairing between fermions: there is no way to 
form a SU(iV) singlet with only two fermions. The only super- 
fluid instability that can be stabilized is a molecular one where 
A fermions form a SU(A) singlet: Mj = c\ ^cj ^ • • • c' N ^ In 
this letter, we shall show, by means of a combination of an- 
alytical and numerical results obtained by the density-matrix 
renormalization group (DMRG) technique lfl5ll . that this MS 
phase emerges in the phase diagram of model ([T]i for U < 
and at at small enough density n. The latter phase is not an 
artifact of the extended SU(AT) symmetry of model (Q~|) — it is 
robust to symmetry breaking terms toward more realistic situ- 
ations. In this respect, we believe that the Hubbard model (fl} 
captures the main generic features responsible for the forma- 
tion of the MS phase. 

Low-energy approach. The low-energy effective field the- 
ory corresponding to the SU(iV) Hubbard chain (QJ can be 
derived, as usual, from the linearization at the two Fermi 
points (±&f) of the dispersion relation of free A-component 
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fermions 11161 Il7ll . The derivation of the low-energy Hamil- 
tonian is straightforward (see for instance Ref. [l8| for details) 
and, away from half-filling, it separates into two commuting 
charge and spin pieces: H = H c + Ti s . This is the famous 
spin-cha rge sep aration which is the hallmark of ID electronic 
systems Ill6lll7ll . Within this low-energy description, the U(l) 
charge excitations are described by a free massless bosonic 
field $ with Hamiltonian density: 

1 



K 



(d x $) 2 + K (d x G) 2 



(2) 



where v and K are respectively the charge velocity 
and the Luttinger parameter. A perturbative estimate 
gives v = v F [1 + U(N - 1)/ (nv F )} 1/2 and K = 
[1 + U(N — 1)/(ttVf)] with the Fermi velocity vp = 
2tasm(kFa) and lattice spacing a. For generic fillings, no 
umklapp terms appear and the charge degrees of freedom dis- 
play metallic properties in the Luttinger liquid universality 
class |fl6l[l7ll . 

The hyperfine spin sector is described by the SU(AT) 
Thirring model which is an integrable field theory 11911 . For at- 
tractive interaction (U < 0) a spectral gap m opens. The low- 
energy spectrum in the hyperfine spin sector consists of N — 1 
branches with masses m r = m sm(irr/N) (r = 1, . . . , N— 1) 
lfl9ll . The dominant instability which governs the physics 
of this phase is the one with the slowest decaying correla- 
tions at zero temperature. Both the 1 -particle Green func- 
tion G(x) = (c' a jCa^+rr) and (onsite) pairing correlations 

P(x) = (cl i c' /3 , i Cf3.i +x c at i +x ) are short-range. On the con- 
trary, the equal-time density correlation N(x) = (riirii+x) 
associated with a charge-density wave (CDW) and the equal- 
time MS correlations M(x) = (MjM} + „) have the following 
power-law decay at long distance 12011 : 

cos(2k F x)x~ 2K/N 

X -N/(2K) for7Veven 



N(x) 
M(x) 

M(x) ~ sm(k F x)x^ K+N2 / K ^^ for N odd. 



(3) 
(4) 
(5) 



We thus see that CDW and MS instabilities compete and 
the key point of the analysis is the one which dominates. At 
issue is the value of the Luttinger parameter K. In particular 
a dominant MS instability requires K > N/2 (K > A/V3) 
for N even (odd, respectively) and thus a fairly large value 
of K which, with only short range interaction, is not guaran- 
teed. However a simple argument suggests that this may be 
realized at sufficiently small density at large negative U. In- 
deed, when n <C 1 and \U\/t ^S> 1, a dilute gas of strongly 
bound A-fermion objects forms and (Q~|) behaves as essen- 
tially free hardcore bosons (A even) or free fermions (A odd) 
with an effective hopping t /\U\ . One can therefore es- 
timate M(x) in this limit as the free bosonic Green function, 
M(x) ~ x^ 1 ! 2 1 17, 2lll . when A is even and, as the free 
fermion Green function, M(x) ~ sm(k F x) / x, when A is 
odd. By comparing with Eqs. (01 |5kwe deduce an upper 
bound for K which is K max = TV 12211 . From the pertur- 
bative estimate we see that K > 1 and K increases with 
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FIG. 1: (Color online) SU(3) model: triplet and density correlations 
vs distance obtained by DMRG with L = 153, n = 1/3 and U/t = 
—4. (a) dominant triplet over CDW correlations can both be fitted 
with K — 2.7. We also see the k F (respectively 2fcp) oscillations of 
M{x) (respectively N(x)). (b) 1-particle Green function G(x) and 
pairing correlations P(x) vs distance. Both are short-range and with 
the same correlation length £ = 0.68. 



\U\, so that there is room to stabilize an MS phase for suf- 
ficiently strong attractive interaction and small density. In ad- 
dition, at zero density, the A^-component Fermi gas with an 
SU(A) symmetry is known to be exactly-solvable by means 
of the Bethe-ansatz approach and bound states of N fermions 
are formed for attractive interaction f23 : <]. Outside these cases 
of infinite attractive interaction or vanishing density, the exis- 
tence and stability of this MS phase stem from the full non- 
perturbative behavior of the Luttinger parameter ifasa func- 
tion of the density n and the interaction U. We shall now 
evaluate numerically this parameter in the simplest odd and 
even cases N = 3, 4 by computing dominant correlations with 
the DMRG technique to conclude on the extension of the MS 
phase. 

Numerical Results. We have performed extensive DMRG 
calculations for both the N = 3 and N = 4 cases and for 
a wide range of densities n and interactions U 124 12511 . We 



show in Fig. |TJa-b) and Fig. |2a-b) our data for N = 3 and 
N = 4 respectively at typical values of n and U = —At. In 
both cases, and in agreement with the low-energy approach, 
we find that a gap opens in the hyperfine spin sector and 
that the one-particle and pairing correlations are always short 
ranged. From Fig. |TJb) and Fig. EJb), one can compute the 
one- and two-particle correlation lengths £ that is expected to 
vary as the inverse gap : £ ~ 1/m. We find that the ratio 
R = m\/m.2 is close to 1 and l/y/2 respectively for N = 3 
and N — 4 as expected from the low-energy approach. In 
contrast, we see in Fig.QJa) and Fig. [2 a) that the density and 
MS correlations N(x) and M (x) display power-law behav- 
ior. Clearly, triplet and quartet correlations dominate over 
CDW at these densities (n = 1/3 for N = 3 and n = 0.5 
for N = 4). The phase diagrams for both SU(3) and SU(4) 
models are presented in Fig. [3] and Fig. [4] which give a map 
of K vs interaction and density. The values of K were ob- 
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tions vs distance obtained by DMRG with L = 128 and U/t — — 4 
at filling n = 0.25. The same K ~ 2.7 is used to give a rough esti- L//t 
mate of the exponent, (b) 1- and 2-particle correlations vs distance. 

Both are short-range and the ratio of the two correlation lengths is FIG. 4: Same as Fig.[3]for the SU(4) case. 



= 0.68 ~ 1/^2. 
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FIG. 3: SU(3) model: phase diagram showing the Luttinger parame- 
ter K vs filling n and interaction U. The grey area is the superfluid 
triplet phase. Lines are guide for the eyes which behavior satisfies 
the perturbative limit close to the point (U/t = 0, n = 0). We also 
plot the perturbative estimate separating both regions (see text), that 
agrees with our numerical finding at small \ U\ but deviates for larger 
values. 



tained from the power-law behavior of the molecular correla- 
tion M(x) using Eqs. (01 [5]). We find that triplet and quartet 
superfluid phases emerge in a wide portion of the phase dia- 
grams (grey area) separated from a CDW phase by a cross- 
over line n c (U). As the density decreases from half-filling 
(n = N/2) to 0, K increases from N/4 to N for moderate 
or strong |£7|/t, as expected from the strong-coupling argu- 
ment. Interestingly enough, the MS phase extends to small 
values of U at sufficiently small densities. We also observe 
that the curve n c (U) is likely to saturate in the strong cou- 



pling limit l25ll . For example in the SU(4) case, the K(n) 
function is almost [/-independent for \U\/t > 2 so that lines 
of equal K are parallel to the U axis. 

Effect of symmetry breaking perturbations. At this point, 
the natural question is whether the molecular superfluid 
phases survives to the breaking of the SU(A r ) symmetry. This 
is an important question since in most of the realistic situa- 
tions, the actual symmetry is expected to be much smaller. 
Part of the answer is given in ID systems by the accepted 
view that, at sufficiently low energies and for generic in- 
teractions, the dynamical symmetry is most likely to be en- 
larged [26]: though the S\J(N) symmetry is not an exact sym- 
metry, it is physically meaningful as an effective low-energy 
theory. As an example, we consider the SU(4) case rele- 
vant for spin-3/2 cold atoms and add to the Hubbard Hamil- 
tonian (Q]) a singlet-pairing coupling Vj^i Pqo i^oo,i where 

P 00,i = 4/2,i C -3/2,i ~ 4/2,/- 1/2,* ■ AS Sh ° wn in Ref.EZ, 

the pairing term reduces the SU(4) symmetry down to SO(5). 
We show typical data for U/t = -4 and V/t = -2 at the 
density n = 1/2 in Fig. [5] We clearly see that the equal-time 
pairing correlation function P(x) = (Pj jPoo.i+a;) admits an 
exponential decay i.e. there is no BCS instability. In contrast, 
quartet correlations are (quasi) long ranged and dominate over 
CDW ones. Remarkably, we observe from Fig. [5]that the gap 
ratio R is very similar to the one for the full SU(4) symmetric 
model when V = 0. This means that the SU(4) model (HJ is 
a very good starting point to explore the main features of the 
quartet phase. Of course, for large negative V, a BCS phase 
does appear 112 8ll but the main point here is to show that the 
quartet molecular phase is not an artifact of the SU(4) sym- 
metry and does exist in more realistic models I29I1 . A more 
detailed study will be presented elsewhere. 

Concluding remarks. We have shown that a quasi-long- 
range general MS phase can emerge in ID for attractive in- 
teractions at low density. This ID phase, characterized by 
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FIG. 5: SO(5) model : (a) Quartet and density correlations for U/t = 
—4 and V/t = —2 at the density n=0.5. The value of the Luttinger 
parameter is K = 2.7. (b) 1- and 2-particle correlations vs distance. 
Both are short-range and the ratio of the two correlation lengths is 
6/6 = 0.72 ~ 1/V2. 

a bound-state made of TV fermions, can be viewed as a ne- 
matic Luttinger liquid and a simple paradigmatic model to 
describe its main physical properties is the attractive SU(7V) 
Hubbard chain (Q]). The triplet and quartet phases in the sim- 
plest N = 3, 4 cases might be explored experimentally in the 
context of spinor ultracold fermionic atoms. As a first step, we 
have assumed here a homogeneous optical lattice and neglect 
in the first approximation the parabolic confining potential of 
the atomic trap. We expect that this potential will not affect 
the properties of this molecular phase at low density. Such an 
effect could be investigated by DMRG calculations for quan- 
titative comparisons [30]. In the context of cold atoms exper- 
iments, the triplet and quartet phases can be probed by radio- 
frequency spectroscopy to measure the excitation gaps of the 
successive triplet-quartet dissociation process. We hope that 
future experiments in ultracold fermionic atoms will reveal the 
existence of these triplet and quartet phases. 
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